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1.0  Introduction 


The  Global  Positioning  System  (GPS)  has  been  used  successfully  for  high  precision 
relative  positioning  for  nearly  a  decade.1  The  subcentimeter  level  accuracy  achieved  is 
attributed  to  the  use  of  reconstructed  carrier  phase  measurements  of  the  signals  received  from 
GPS  satellites.  Additional  improvements  in  accuracy  have  been  obtained  through  the  use  of  the 
Precise  (P)  code  pseudorange  information.2-3  This  improvement  is  significant  if  the 
observation  period  is  short,  making  ambiguity  resolution,  without  the  P  code  delay 
observations,  impossible.  The  lower  accuracy  Coarse/Acquisition  (C/A)  code  has  not  been 
widely  used  together  with  phase  observations.  The  purpose  of  this  study  is  to  evaluate  the 
possible  improvements  in  relative  positioning  which  might  result  from  addition  of  C7A  code 
measurements  to  dual  frequency  phase  observables.  This  work  is  motivated  by  changes  which 
will  take  place  when  the  new  constellation  of  Block  II  satellites  become  active.  All  present 
work  has  been  done  using  the  aging  constellation  of  Block  I  satellites.  These  satellites  were 
designed  to  verify  the  accuracy  of  GPS  positioning  for  defense  applications,  and  to  support 
development  of  GPS  user  equipment.  This  constellation  consists  of  six  operational  satellites, 
in  two  orbital  planes,  oriented  to  provide  maximum  coverage  over  North  America.  In  the 
Block  I  satellites  none  of  the  codes  are  encrypted.  As  the  GPS  system  becomes  operational 
with  launching  of  the  Block  II  satellites,  the  P  code  will  be  encrypted  for  national  security 
reasons.  The  P  code  provides  the  highest  accuracy,  and  unauthorized  users  possessing  the 
code  could  use  the  knowledge  to  deceive  authorized  users.  The  introduction  of  advanced  GPS 
receivers  provides  a  second  motivation  for  this  study. 


'Counselman  III,  C.  C.,  ct.  af  .  Accuracy  of  Baseline  Delerminanons  bv  MITES  Assessed  by  Comparison  wuh 
Tape.  Theodolite,  and  Geodimeter  Measurements EOS.  the  Transactions  of  the  American  Geophysical  i'nton.  sol  62. 
p  260.  April  28.  1981. 

2Wubbcna.  G  .  A  Schuchardt,  and  G  Seebcr,  Multistation  positioning  results  with  TI4100  GPS  receivers  ;n 
geodetic  control  networks  '.  Proc  Fourth  International  Geodetic  Symposium  on  Satellite  Positioning,  soi  2.  pp  9*1. 
978.  University  of  Texas.  Austin,  1986 

3Lichten,  S  M  ,  and  J  S  Border,  Strategies  for  High  Precision  Global  Positioning  System  Orbit  Determination  , 
Journal  of  Geophysical  Research,  vol  92.  pp  12751  12762.  1987 


Geodetic  GPS  receivers  with  the  capability  to  measure  precise  dual  frequency  carrier 
phase  as  well  as  C/A  code  delay  have  been  introduced.  The  Mini- Mac™  was  the  first  receiver 
to  enter  the  market.  Several  receivers  with  similar  capabilities  are  also  being  marketed  or  have 


been  announced. 

The  advantage  of  utilizing  code  delay  information  is  that  it  is  an  unambiguous, 
relatively  unbiased  measurement  of  the  satellite-receiver  range.  The  earner  phase  measurement 
is  more  precise,  but  contains  a  bias  which  must  be  determined  to  recover  range  informauon. 
The  drawback  of  using  delay,  or  pseudorange,  information  is  the  relatively  large  amount  of 
noise  present  in  it.  A  further  disadvantage  in  using  the  C/A  code  is  its  availability  at  only  a 
single  frequency,  making  removal  of  the  ionospheric  delay  difficult  It  is  hoped  that  by 
minimizing  ionospheric  contnbution  and  noise,  C/A  code  can  enhance  phase  observations. 


Mini-Mac  is  a  registered  trademark  of  Aero  Sersicc  Division. 


Western  (ieophysical  Company  of  America 


2.0  Background  material 

This  section  will  first  describe  the  two  types  of  GPS  observables  which  are  used  for 
positioning.  These  observables  are  the  phase  and  pscudorange  measurements.  Next,  the 
errors  which  corrupt  these  measurements  will  be  discussed. 

2.1  GPS  Signal  Structure 

The  GPS  satellites  transmit  positioning  signals  in  two  radio  frequency  bands,  the  L 1 
band  and  the  L2  band,  with  carrier  frequencies  of  1575.42  MHz  and  1227.60  MHz 
respectively.  Both  periodic  carrier  waves  are  phase  modulated  by  pseudorandom  codes.  The 
codes  appear  as  random  noise  to  those  not  possessing  the  keys.  Two  different  pseudorandom 
codes  are  used  by  the  GPS  system.  A  so-called  Coarse/Acquisition  (C/A)  code  modulates  the 
LI  carrier.  A  second  code,  known  at  the  Precise  (P)  code,  modulates  both  earner  waves.  As 
their  names  implies,  the  P  code  can  provide  gi eater  positioning  accuracy  than  the  C/A  code. 
The  C/A  code  is  used  to  acquire,  or  lock  on  to  the  GPS  signals.1  The  modulation  suppresses 
the  earner  wave;  however  a  earner  wave  can  be  reconstructed.  The  wavelength  of  the  earner 
waves  are  19  cm  for  LI  and  24  cm  for  L2.  Range  information  can  be  derived  by 
measurements  of  the  reconstructed  carrier  phase  and  the  code  modulation  delay. 

2.2  Phase  observable 

The  advantage  of  measuring  the  carrier  phase  is  the  precision  with  w hich  the 
measurement  can  be  made.  Using  a  well  designed  antenna  and  receiver,  the  phase  can  be 
measured  to  one  percent  of  a  cycle,  which  corresponds  to  about  2  mm  of  range. 

High  accuracy  relauve  positioning  is  readily  obtained,  using  differenced  phase 
measurements.  In  relative  positioning  the  position  of  one  station  is  know  n,  and  the  position  of 
the  second  station  is  determined  with  respect  to  the  fixed  station.  The  position  solutions  are 

'Spilker,  J.  J  .  CPS  signal  structure  and  performance  characteristics  '.  (ilobal  Positioning  Ssstem,  The  Institute 
of  Navigation,  vol  1,  pp.  29  54,  1980 
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generally  not  available  in  real  time.  Typically,  observations  taken  over  several  hours  are 
processed  to  determine  baselines,  with  phase  measurements  typically  taken  every  six  minutes. 

2.2.1  Codeless  measurement 

Several  methods  are  available  to  measure  the  phase  of  the  reconstructed  carrier.  If  the 
receiver  has  knowledge  of  the  modulating  codes  used  by  the  GPS  satellite,  the  carrier  wave 
can  be  reconstructed.  In  this  case  the  wavelengths  of  carrier  are  19  cm  and  24  cm.  for  the  L 1 
and  L2  bands  respectively.  This  is  the  method  used  m  the  TI  4100  receiver  in  both  LI  and  L2 
bands.  If  the  receiver  does  not  possess  knowledge  of  the  codes,  the  second  harmonic  of  the 
incoming  earner  is  reconstructed.  The  codeless  technique  results  in  measurements  made  with  a 
reconstructed  carrier,  having  wavelengths  of  9.5  cm  and  12  cm  respectively  for  the  LI  and  L2 
bands.  The  shorter  carrier  wave  observations  are  obtained  from  the  Macrometer  11  ^  at  both 
LI  and  L2  frequencies.1  The  Mini-Mac  receiver  is  designed  to  use  LI  code  and  no  L2  code.2 

2.2.2  Integer  ambiguity 

The  one-way  phase  observable  is  a  precise  measurement  of  the  fractional  pan  of 
reconstructed  carrier  wave.  Since  this  measurement  is  corrupted  by  errors,  the  unknown 
ambiguity  is  not  an  integer.  Bv  differencing  observations  between  satellites  and  receivers,  the 
errors  are  significantly  reduced.  The  ambiguity  parameter  of  the  differenced  observable  is  an 
integer.  By  determining  the  ambiguity  parameter,  uncertainties  in  positioning  can  be  reduced. 
The  process  of  resolving  the  ambiguity  is  known  as  bias  fixing 


TV 

Micrometer  LI  is  a  registered  trademark  of  Aero  Service  Division,  "■'tsitm  Geophysical  Company  of  America 

'King.  :>  .W  ,  E  G.  Masters.  C  Ri/os.  A  Siolr.  and  i  Collins,  Surveying  wiih  ( ,PS.  Monograph  No  9.  School 
Surveying.  University  of  New  South  W  ales.  Kensington.  NSW,  Australia,  November  19X5 

2Ladd,  J  W  .  R  G.  Welshe,  .Mini  Mac'*  A  New  Generation  Dual  Band  Surveyor  Prut  Fourth  InlernaUuruii 
Geodetic  Symposium  on  Satellite  PtiMiu-ntng.  vul  <,  pp  475  4X7.  L'mversiiy  of  Texas.  Austin.  19X6 


4 


2.3  Pseudorange  observable 

The  second  type  of  position  information  available  is  the  time  delay  of  the  code 
modulation.  The  code  observables  serve  to  determine  the  unambiguous  range  of  the  satellites 
and  to  identify  each  satellite. 

The  satellite-receiver  range  is  determined  by  measuring  the  time  of  the  reception  of  the 
code  modulated  signal.  The  satellites  transmit  unique  codes  which  are  correlated  with  a  replica, 
generated  simultaneously  in  the  receiver.  The  receiver  code  is  delayed  to  maximize  the 
correlation  with  the  received  signal.  The  added  delay  is  a  measure  of  the  satellite-receiver 
range.  This  measured  range  is  termed  pseudorange,  since  it  is  corrupted  by  clock  errors.  If 
receiver  and  satellite  clocks  are  synchronized,  three  satellites  can  accurately  determine  the 
receiver  position  in  three  dimensions.  If  receiver  and  satellite  clocks  are  not  synchronized,  four 
satellites  are  needed.  The  redundant  information  is  used  to  solve  for  relative  receiver-satellite 
clock  error. 

The  coded  modulation  makes  them  resistant  to  intentional  and  unintentional  jamming. 
To  prevent  all  the  satellites  from  interfering  wit)  each  other,  the  codes  chosen  have  the  property 
of  being  orthogonal  to  each  other  Any  cross  correlation  between  any  two  codes  will  be  zero, 
and  only  an  autocorrelation  results  in  a  nonzero  value.  Different  sections  of  the  same  code  are 
also  orthogonal  to  each  other.  This  property  is  helpful  against  multipath  propagation  errors,  as 
will  be  discussed  in  the  section  on  errors.1 
2.3.1  Coarse/Acquisition  (C/A)  Code 

The  C/A  code  occupies  a  bandwidth  of  approximately  1  MHz.  The  axle  consists  of 
approximately  KP  binary  bits  which  is  modulated  at  a  rate  of  1  Mbps,  reuniting  in  a  duration  of 
1  msec  per  cycle  of  the  C/A  axle.  This  code  is  used  for  acquiring  the  GPS  signal  and 
determining  the  satellite-receiver  propagation  delay.  Information  concerning  clock  corr-ction, 
satellite  ephemens,  and  system  messages  also  modulates  the  Cl  earner  wave.1 

1  Spilker.  J  J.  Ci PS  signal  s (nature  ar.«l  performance  characteristics  < »7- P  s.U.-nint  .Vv«.v*n,  The  1:>;.:.::j 
•jf  Navigation,  vol  1 .  pp  29  Vl. 
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The  C/A  code  repeats  itself  evw.y  millisecond,  implying  that  the  receiver  and  satellite 
clocks  must  be  synchronized  to  within  a  millisecond,  or  the  a  prion  receiver  position  must  be 
known  to  within  300  km  to  determine  unambiguous  range  measurements.  A  typical  C/A  code 
only  receiver  can  determine  posinon  to  within  30  m. 

2.3.2  Precise  (Pi  code 

Unlike  the  C/A  code,  which  will  be  available  to  everyone,  the  P  code  will  be  restricted 
to  national  security  applicanons,  when  the  Block  II  satelli  *s  are  operational.  The  P  code 
access  is  restricted  to  prevent  the  jamming  of  the  GPS  signals.  A  ground  or  satellite  transminer 
broadcasting  the  P  code,  near  a  receiver,  will  make  accurate  posmoning  impossible.  The  P 
code  occupies  a  bandwidth  of  approximately  10  MHz.  The  total  code  length  is  267  days,  with 
each  satellite  assigned  a  unique  one  week  segment  of  the  code.  Unlike  the  C/A  code,  the  P 
code  is  broadcast  on  both  the  LI  and  L2  bands,  enabling  the  ionospheric  effect  to  be  reduced, 
and  permitting  higher  accuracy.  With  the  P  code,  three  dimensional  position  can  be  determined 
to  within  10  meters  in  real-time.1 

2.4  Observation  errors 

Several  factors  contribute  to  degradation  of  both  phase  and  pseudorange  measurements, 
but  to  different  degrees.  These  contributions  are  discussed  in  the  following  sections. 

2.4.1  Oscillator  phases 

The  largest  source  or  error  on  one-way  observables  are  the  local  oscillator  phases, 
present  in  each  satellite  and  receiver.  These  errors  can  be  readily  removed  by  differencing 
observations  between  receivers  and  satellites.  Uach  of  t he  GPS  satellite  is  equipped  with 
Cesium  oscillators,  where  as  the  GPS  receivers  are  equipped  with  relatively  inexpensive,  and 
less  accurate,  quartz  oscillators.  Both  receiver  and  satellite  oscillator  phases  will  not  permit  the 
integer  ambiguity  to  be  determmed,  unless  these  errors  are  significantly  reduced. 

'  Spilker,  J  J,  GPS  signal  strut  tu  re  anU  pet  formant  e  ,  ara,  icr  ..its  .  On -ttj  ISsiuonmg  System.  The 
Institute  of  Navigation  vol  1,  pp  54.  I'rsO 
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2.4.2  Satellite  Ephemeris 

The  inability  to  accurately  predict  the  orbit  of  the  GPS  satellites,  also  limits  the 
positioning  accuracy.  These  error  are  also  readily  reduced  by  forming  differences.  The 
satellites  are  a  finite  distance  form  the  receivers,  and  therefore  the  satellite  ephemens  errors  are 
passed  onto  baseline  measurements.  The  difficulties  in  orbit  prediction  arise  from  the  inability 
in  accurately  modeling  all  the  forces  acting  on  the  satellites,  and  being  unable  to  precisely 
measure  the  initial  conditions  of  the  satellites.  Ihe  ability  to  measure  baselines,  with  fixed 
orbits,  depends  on  the  following  empirical  relationship:1 

baseline  error  _  orbit  error 
baseline  length  height  of  satellite 

With  orbit  errors  of  tens  of  meters  and  the  height  of  GPS  satellites  at  20,183  km,  the  accuracy 
achievable  is  about  one  part  in  107.  To  improve  baseline  accuracy,  knowledge  of  the  satellite 
orbits  must  be  improved.  It  has  been  shown  that  the  ability  to  fix  biases  on  smaller  baselines 
enhances  the  ability  to  determine  orbits.  Fixing  biases  on  a  baseline  reduces  the  position 
uncertainty  of  the  receivers,  permitung  higher  accuracy  orbit  determination.2 

2.4.3  Propagation  delay 

Significant  error  in  range  measurement  is  introduced  by  the  environment, 
through  which  the  GPS  signals  travel.  One  environmental  error  source  is  the  ionospheric 
delay,  whose  effect  is  inversely  proportional  to  the  square  of  the  frequency  for  code 
observables,  and  inversely  proportional  to  frequency  for  phase  observables.  This  property 
enables  the  the  ionospheric  contribution  to  be  removed  by  making  measurements  at  two 
frequencies  and  forming  a  linear  combinauon  of  two  measurements.  This  correction  wall 
discussed  in  more  detail  in  Section  4. 1  The  ionosphere  has  different  effects  on  the  phase  and 

•Colombo.  OL.  fphcmcrts  terrors  of  Cl'S  Satellites  .  Bulletin  (leodeuque .  soi  60,  pp  64  84.  WS6. 

‘Abbot.  RI.CC  Cnunselman  111  or  1  S  A  Courevik h.  Cl'S  Orbit  Determination  Bootstrapping  to  Rosotse 
f  airier  Pha.se  Ambiguity  .  presented  at  the  hy.H  /nierncit.iuwl  Ctthlel i<:  .S' vmposiu/n  on  Sau-U.'.e  /V.vir.,  nmg,  Sow. 
Mmmi  State  L'mvrrs.tv.  Lii  Cruses.  March  1<  I7.  Idg<J 
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pseudorange  observables.  The  ionosphere  advances  the  phase  as  the  signal  propagates  through 
it,  thereby  decreasing  the  measured  range  to  the  satellite.  The  ionospheric  delay  for  the 
pseudorange  observable  is  opposite  in  sign  to  delay  for  the  phase  observable.  The 
pseudorange  observable  is  delayed  by  the  ionosphere,  resulting  in  larger  range  measurement. 
Ionospheric  delay  can  range  from  several  meters  during  the  night  to  tens  of  meters  during  the 
day.  The  ionospheric  delay  is  also  dependent  on  solar  activity,  satellite  elevation,  and  latitude 
of  the  GPS  receiver.  The  ionospheric  delay  is  large  during  periods  of  strong  solar  activity,  low 
satellite  elevations,  and  low  latitudes.  Another  source  of  propagation  error  is  the  troposphere. 
Both  water  vapor  and  other  dry  atmospheric  components  further  delay  the  GPS  signals. 
Tropospheric  delay  can  be  readily  modeled  by  using  basic  atmospheric  measurements  at  the 
surface,  such  as  temperature,  pressure,  and  relative  humidity.  Unmodeled,  the  troposphere 
can  produce  errors  on  the  order  of  tens  of  meters  at  low  satellite  elevations  (  »  '  5°  ).  Using  the 
three  surface  measurements,  the  tropospheric  delay  can  be  reduced  to  approximately  10 
centimeters  or  less.1 
2.4.4  Measurement  errors 

Two  sources  of  error,  which  are  functions  of  receiver  and  antenna  design,  excluding 
clocks,  are  thermal  noise  and  multipath.  Thermal  noise  in  a  well  designed  GPS  receiver  is 
limited  to  1/2  %  -  1  %  of  the  wavelength  of  the  earner  wave  or  the  chipwidth  of  the  codes. 

Multipath  propagation  is  the  reflected  signal  which  reaches  the  antenna  via  a  path  which 
is  not  the  line  of  sight  path.  Reflected  signals  interfere  with  the  direct  signal,  corrupting  the 
measurements.  The  effect  of  multipath  propagation  has  a  different  ef  fect  on  each  observable. 
With  the  use  of  a  metallic,  horizontal  ground  plane  below  the  antenna  dements,  multipath 
effects  on  phase  observations  can  be  reduced  to  subcentimeter  levels.2  For  the  pseudorange 


'Spilker,  J.  J .  GPS  signal  structure  and  performance  characteristics”.  Global  Prmiiorung  System,  The  Institute 
of  Navigation,  vol.  1.  pp.  2934,  ldgo 

“■Courset/nan  III.  C  C  .  S  A  Ciimrev itch.  Miniature  Interferometer  Terminals  for  Earth  Surve>ing:  Ambiguity 
and  Multipath  with  Global  Positioning  System  If.f.P  l ransom-turns  .>n  Geoscience  anj  Remote  Sensing,  vol  GE  19. 
no.  4,  pp,  244-232.  October  19H! 
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observable,  the  range  error  due  to  multipath  propagation  is  limited  to  the  chip  width  of  the  code 
being  received.  If  the  difference  between  the  direct  and  reflected  signal  path  is  greater  than  one 
chip  width,  the  receiver  automatically  rejects  this  signal  due  to  the  orthogonality  property  of  the 
codes.  On  the  other  extreme  is  the  case  when  the  reflected  signal  path  is  approximately  equal  to 
the  line  of  sight  path.  In  this  case  the  multipath  effect  is  on  the  order  30  cm  for  F  code.2 
Multipath  propagation  error  is  expected  to  be  on  the  order  of  1  m  for  one-way  OA  code 
observations.  This  is  the  case  for  geodetic  receivers,  where  the  antenna  is  mounted  1-2  m 
above  the  ground  When  the  antenna  is  this  close  to  the  ground  and  stationary,  multipath 
propagation  errors  exhibit  quasipenodic  variations,  due  to  the  satellite's  motion  across  the  sky. 
The  period  of  these  variations  is  on  the  order  of  1 5  minutes.  This  periodicity  implies  that  a 
partial  averaging  of  these  errors  occur  over  time.1 


'Couaselmin  III,  C  C  .  personal  communication  (  MIT). 
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3.0  Differencing  observables 

Oscillator  phase,  ionospheric,  and  tropospheric  contributions  must  be  significantly 
reduced  to  permit  the  integer  ambiguity  to  be  determined.  These  effects  can  be  reduced  by 
forming  differences  among  a  common  set  of  observations,  canceling  common  mode  errors 
present  in  the  observables. 

3.1  Single  Differencing 

If  two  receivers  simultaneously  observe  the  same  satellites,  the  observations  from  both 
receivers  can  be  differenced.  The  new  observable  is  known  as  a  between  receiver  single 
difference  observable.  The  geometry  for  single  differencing  is  shown  in  Figure  3.1. 


Figure  3. 1  Satellite-receiver  configuration  for  a  single  difference  observable 


In  performing  the  single  difference,  the  measured  phase  observations  between  the  satellite  and 
each  receiver  is  differenced.  In  the  process  of  being  differenced,  errors  common  to  both  range 
measurements  cancel  each  other.  Phis  can  be  seen  by  forming  a  simplified  theoretical  equation 
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for  a  one-way  observable.  The  phase  that  is  measured  at  the  receiver  is  a  function  of  several 
factors  shown  in  Equation  3. 1.  The  notation  used  is  from  Footnote  1  given  below. 


0k  =  0  +  ok  - 


rk  +  nk 


(3.1) 


i  -  satellite  index 

k  -  receiver  index 

<5  k'  -  measured  phase  of  signal  at  receiver  k  from  satellite  i 
0  1  -  satellite  i  oscillator  phase 

0  k  -  receiver  k  oscillator  phase 
k  -  wavelength  of  the  reconstructed  carrier  signal 

r  k‘  -  actual  range  between  receiver  k  and  satellite  i 

n  k'  -  integer  number  of  cycles  between  receiver  k  and  satellite  i 
at  the  initial  epoch 
e  k‘  -  errors  in  measured  range 

-  (i.e.  ionosphere,  troposphere,  and  satellite  ephemens) 


The  phase  of  the  carrier  signal  increases  as  the  range  decreases,  and  decreases  for  increasing 
range,  resulting  in  a  negative  sign  in  front  of  the  range  term.  Using  an  additional  one-way 
observable  by  a  second  receiver  (q),  a  single  difference  observable  can  be  formed: 


q  -  receiver  index 


Equation  3.2  can  be  simplified  by  denoting  a  single  difference  by  (A): 


*B«utler,  G  ,  et.al  ,  l  sing  the  Global  Positioning  System  ;GPSl  for  High  Precision  Geodetic  Suness 
Highlighu  and  Problem  Areas  ’,  IFFF.  P LA  VS'  M6  Positions  LixMion  j nJ  Symposium,  pp  24?  250. 

November  1986. 


1  1 


(3.3) 


A<t>*q=!*k-*q)  ArU  +AnU+AEicq 

One  of  the  primary  benefits  of  single  differencing  is  that  satellite  oscillator  phases  cancel.  If  the 
receivers  are  not  widely  separated,  both  the  receivers  essentially  "see"  the  same  ionosphere  and 
troposphere;  therefore,  these  errors  are  significantly  reduced.  The  satellite  ephemeris  errors  are 
also  reduced.  Reductions  of  the  ionospheric  contribution  and  the  ephemeris  errors  are  a 
function  of  baseline  length. 

3,2  Double  Differencing 

The  procedure  of  differencing  the  observables  can  be  expanded  further,  to  cancel 
receiver  oscillator  phases.  In  the  single  difference  observable,  satellite  oscillator  phases  cancel, 
but  the  difference  of  the  two  receiver  oscillator  phases  remain.  By  adding  a  satellite  to  the 
scenario  shown  in  Figure  3.1,  a  additional  single  difference  observable  can  be  formed.  Both 
these  differences  do  not  contain  satellite  oscillator  phases,  but  they  do  contain  identical  receiver 
oscillator  phases.  If  one  single  difference,  formed  with  satellite  i,  is  differenced  with  another 
single  difference,  formed  with  satellite  j,  the  resulting  equation  is  shown  below: 


j  •  satellite  index 

The  receiver  oscillator  phases  now  cancel  each  other  completely.  The  resulting  observable  is 
known  as  a  double  difference  observ  able,  and  it  contains  no  satellite  or  receiver  oscillator 
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phases.  If  (AA)  is  used  to  d  note  a  double  difference  quantity,  Equation  3.4  can  be  simplified 
to: 


AAOkJq  =  -  -  A  A  r  kJq  +  A  A  nkJq  +  AAek 


1  j 


(3.5) 


This  configuration  for  forming  a  double  difference  observable  is  shown  in  Figure  3.2. 
Because  the  oscillator  phases  have  been  canceled,  the  double  difference  integer  ambiguity 
parameter  can  be  estimated  as  an  integer.  In  practice,  double  difference  observables  are 
processed  to  estimate  ne  integer  ambiguity,  and  to  determine  the  baselines  between  receivers. 


Figure  3.2  Satellite-receiver  configuration  for  a  double  difference  observable 

Although  the  double  differencing  scheme  is  straightforward  for  two  receivers  and  two 
satellites,  it  becomes  more  complicated  when  more  than  two  satellites  are  involved.  Earlier 
work  involved  picking  one  satellite  as  the  reference  satellite.  This  satellite  was  usually  the  first 
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satellite  visible  in  the  observing  period.  Observations  from  the  remaining  satellites  would  be 
differenced  with  the  reference  satellite.  The  drawback  to  this  procedure  was  that  the  reference 
satellite  appeared  in  the  double  differences  more  than  the  other  satellites,  and  was  given  more 
weight  in  processing.  More  information  was  extracted  from  the  reference  satellite,  than  from 
the  other  satellites.  With  all  the  satellites  having  the  same  quality  of  orbits,  it  was  incorrect  to 
give  more  weight  to  one  particular  satellite.1  This  procedure  for  double  differencing  also 
results  in  observables,  that  no  longer  have  independent  errors.  In  forming  linear  combinations 
(i.e.  forming  differences)  of  the  observables,  the  observable  errors  become  correlated. 

The  solution  to  this  problem  of  improper  weighting,  and  to  create  uncorrelated 
observables  is  to  orthogonalize  the  set  of  single  differences  formed  at  each  epoch.  This  can  be 
accomplished  using  Gram-Schmidt  orthogonalization.2  The  general  equation  which  forms 
equally  weighted,  independent  double  difference  observables  given  a  set  of  single  differences 
at  each  epoch,  using  the  notation  previously  developed  is:3 


IAA  <J>k  q  = 


i  ♦  l 


1 

1 


m  =  1 


(3.6) 


IAA  -  independent  double  difference  observable 
i  -  loops  through  1  to  (n  -  1 )  satellites  at  each  epoch 


'Abbot.  R.I.,  personal  communication  (  M  I  T  V 

2Strang.  G.,  Linear  Algebra  and  ilt  Applications,  Academic  Press,  Inc  .  New  York.  NY.  pp  128  131.  1980. 

3Remondi.  B.W  ,  "Using  the  Global  Positioning  System  (GPS)  phase  observable  for  relative  geodesy:  modeling, 
processing,  and  results'",  PhD  Dissertation,  Center  for  Space  Research,  The  University  of  Texas  at  Ausun,  1984 
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4.0  Processing  algorithm 

The  methods  by  which  errors  in  the  observables  are  reduced  are  discussed  in  Section 
4.1.  Later  sections  develop  the  normal  equations,  through  which  the  observations  are 
processed.  Normal  equations  for  LI  and  L2  phase  observables  will  be  initially  developed,  and 
the  same  technique  will  be  expanded  to  include  the  CJA  code  observable.  Although  the  double 
difference  notation  (AA)  is  not  used  in  this  section,  all  the  observations  and  parameters  are 
double  differenced  before  being  processed. 

4.1  Processing  LI  and  L2  observables 

The  integer  ambiguity  is  estimated  by  first  predicting  the  satellite-receiver  range,  using  a 
priori  station  coordinates  and  integrated  satellite  ephemerides.  From  this  predicted  range,  the 
expected  phase  observables  are  calculated  and  subtracted  from  the  measured  phases,  to  form 
the  pre-fit  residuals.  If  the  theoretical  model  completely  and  accurately  simulates  reality,  the 
prefit  residuals  will  be  zero.  The  model  can  than  be  used  to  determine  the  integer  ambiguity. 

In  practice,  the  prefit  residuals  are  never  zero.  The  prefit  residuals  are  double  differenced  and 
processed  through  a  maximum-likelihood  estimation  algorithm  to  estimate  the  integer 
ambiguity,  the  ionospheric  contribution,  and  other  parameters.  1  he  maximum  likelihood 
function  is  given  by  Equation  4.1. 


L  (x  , . xn;  a)  =  P|f  (Xi ;  ct) 

1  (4.1) 

f  (  )  -  frequency  function 
n  -  number  of  observables 

\i  -  independent  random  observations 

ci  -  unknown  parameter 


The  function  (f)  is  frequency  function  of  known  random  variable  x,  and  unknown  parameter 
a,  which  can  be  a  function  of  several  additional  parameters.  The  likelihood  function  represents 
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the  joint  probability  of  observing  the  set  of  observations  actually  made.  The  goal  is  to 
determine  a  ,  which  maximize  the  value  of  the  likelihood  function,  given  n  observations  (x,j. 
By  maximizing  the  likelihood  function,  we  a re  maximizing  the  probability  that  values  of  x, 
observed  are  the  most  likely  ones.  If  random  observations  have  normal  distributions,  with 
mean  of  a  and  a  standard  variation  of  Cj,  the  frequency  function  and  the  likelihood  function 
can  be  written  as  shown  below. 


f (X; ;  a)  = 


(4.2) 


,  1  -rlf— ) 

L (xj,...,xn ;  a)  =  -  e  ,  \  I 

n  n 

Urc)2  no, 

i 

a,  -  assumed  a  priori  standard  de\  ution 
for  each  observation 


(4.3) 


The  exponent  must  be  minimized  to  maximize  the  likelihood  function.  This  exponent  is  the 
sum  of  the  squares  of  the  prefit  residual,  or  chi-square  (x2): 


2 


n 


x-Z 

i 


■> 


(4.4) 


The  goal  is  to  determine  which  a  minimizes  X2-  thereby  maximizing  L.1 

If  only  LI  and  L2  phase  observables  are  used,  the  resulting  x2  is  shown  below. 


fSoiiow»y,  C.B  .  Elements  of  the  I  heo 
December  1964. 


r y  of  l)rbtt  Deu'rmiruiiujn.  JP1.  Engineering  Planning  Document.  No 
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2 
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Nepochs 


2 

X, 


(4.5) 


,2  - 

(  J  k  \ 

On  ~  nLl  -  ^geom  ~  *0  1 

-  - - — - +  - 

i  0l2  ~  nL2—  8  0geom  — 

°L!  °L2  (4.6) 


Oli  -  LI  phase  (LI  cycles) 

-  L2  phase  (L2  cycles) 

g  -  ratio  of  L2  and  LI  frequencies 

-  L2/L1  =  60/77 

nu  -  LI  integer  ambiguity 
nL2  -  L2  integer  ambiguity 

dgeom  -  satellite-receiver  range  in  units  of  phase  (LI  cycles) 

-  includes  atmospheric  contribution 

-  function  of  several  parameters 

k  -  theoretical  ionospheric  phase  at  LI  frequency 

-  assumed  independent,  for  every  epoch  and 
satellite-receiver  combination  (LI  cycles) 

Oli  -  assumed  a  priori  standard  deviation 
for  LI  phase  observable  (LI  cycles) 

Ol2  •  assumed  a  priori  standard  deviation 
for  L2  phase  observable  (L2  cycles) 


Using  the  notation  developed  in  Section  3.1,  dgcum  can  be  defined  as: 


A  _ 

4 gcom  ~  ~ 

(4.6a) 


By  minimizing  j}  with  respect  to  the  ionosphere  (k)  and  geometry  parameters  (dgeom),  the 
value  of  the  parameters  are  given  below 
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k  =  _8_  (o^-g^n) 

1  -  8  (4.7) 

^geom  ~  7  (  ^L1  §  C*L2  ) 

1_g‘  (4.8) 

Both  parameters  of  linear  combinations  of  the  LI  and  L2  observables.  The  integer  bias 
parameters  (nu  and  n|_2)  have  been  absorbed  into  its  respective  phase  observable  (i.e.  Q'li  - 
Oli  -  oli  )•  The  expression  for  the  ionosphere  (k)  is  known  as  the  geometry-free  observable, 
since  it  contains  no  geometrical,  or  range,  information.  Equation  4.8  is  known  as  the  "LC" 
observable,  and  is  free  of  ionospheric  effects.  It  is  shown  below. 

‘♦’bC  =  $  Ll  “  2  (  ^  L2_  S  ^  LI  ) 

1-g 

=  0Ll-k  (4.9) 

°ll  =  0Li“nLi  (4.9a ) 

Olc  -  ionosphere- free  phase  observable  (Ll  cycles) 


Equation  4.9  shows  that  the  expression  for  dgcom  is  the  Ll  phase  observable  with  the 
ionosphere  removed. 

Instead  of  using  the  LC  observable,  an  alternative  approach  is  taken  to  permit  the  phase 
cycle  ambiguity  to  remain  an  integer.  This  method  was  developed  by  Dr.  Sergei  Gourevitch  at 
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MIT.  Chi-square  is  rewritten  to  include  a  constraint  on  the  ionospheric  contribution.  Equation 
4.6  can  be  rewritten  as  shown  below. 


dli)n  -  theoretical  ionospheric  contribution  to 
the  LI  carrier  phase  (LI  cycles) 

-  average  value  of  the  ionosphere 

olfin  -  assumed  a  prion  standard  deviation  for 

the  ionospheric  contribution  (LI  cycles) 

-  ionospheric  constraint 

It  is  now  clear  that  Equauon  4  6  is  a  special  case  of  Equation  4.10,  when  al0T1  approaches 
infinity.  On  the  other  extreme  is  the  case  when  c,on  approaches  zero,  implying  that  the 
ionospheric  contribution  is  accurately  known,  and  the  L.1  and  L2  observables  are  processed 
independently  If  no  a  prion  value  for  the  ionosphere  is  available  than  dlon  is  set  to  zero.  A 
ionosphere  model  was  not  used  determine  dum  in  this  study,  and  it  was  assumed  to  be  zero. 
The  ionosphere  (k)  in  Equation  4.10  is  removed  implicitly,  for  every  epoch  and  satellite- 
receiver  combination.  This  is  accomplished  by  minimizing  x~  with  respect  to  k,  and 
substituting  the  expression  for  k  back  into  x In  doing  this,  the  value  of  the  ionosphere  is 
being  estimated  using  LI  and  L2  phase  observables.  The  ionosphere  (k)  is  implicitly 
determined  to  reduce  the  number  of  parameters  that  are  solved  for  in  the  normal  equations. 


Th.  process  of  forming  the  normal  equations  from  Equation  4. 10  will  now  be 
discussed.  The  goal  is  rewrite  '/}  in  a  quadratic  form,  from  which  the  normal  equations  can  be 
readily  derived.  The  quadratic  form  is  shown  below. 


Tr 


Xi  =(  V]  *[  W]  [  V] 


V]  = 


*^geom 

[♦  1 
v  gcom 

^  ion 

0 

Mon  J 

(4.11) 


(4.12) 


Qgeom  -  ionosphere- free  observable 
0,on  ■  geometry-free  observable 


The  weight  matrix  [  V/  )  contains  information  on  the  a  priori  errors  in  the  LI  and  L2  phase 
observables  and  the  ionosphere. 


W  i  , 

W  ,  1 

i  i 
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vv  2  1 

w  2 

(4  13) 
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geom 

_  r  -p  ^ 

^Ll 

0 

ion 

^L2 

(4.14) 


The  transformation  matrix  (  T  |  transforms  1.1  and  L2  observables  to  ionosphere -free  and 
geomt  -y -free  observables 


20 


Finally  x2  can  be  fully  expanded: 


*■*  gcom 
di„„ 
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d  gcom 
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W  2  ] 
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The  terms  in  the  weight  and  the  transformauor.  matrix  were  determined  by  equating  Equation 
4.16  with  Equation  4. 10,  with  the  ionosphere  tern  (k)  implicitly  solved  for.  The  appropriate 
partial  derivatives  were  taken  of  both  sides,  isolating  the  desired  element  in  each  matrix: 
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The  one  half  term  accounts  for  the  factor  two  which  results  from  taking  a  derivative  of  a 
quadratic  term.  Only  the  subscript  1  on  Chi-square  (x:)  represents  the  i-Jl  epoch.  These 
calculations  were  carried  out  by  Dr.  Gourcvuch,  using  Macsyma.  The  full  expressions  for  the 
transformanon  and  weigh;  matrices  are  given  in  Appendix  A.  By  examining  the  transformation 
matrix,  its  dear  that  the  transformed  observables  are  the  same  geometry -free  and  ionosphere- 
free  observables  developed  in  the  previous  example.  The  new  observables  are  a  simple 
rotation  of  the  old  observation  vector.  Since  the  observables  are  processed  using  linearized 
least  squares,  the  geometrical  and  the  ionosphere  models  are  calculated  about  a  a  priori  value. 
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-  a  priori  value  of  satellite  range  (1-1  cycles) 

-  a  priori  value  of  the  ionospheric  contribution  (Ll  cycles) 

-  geometrical  parameters 

-  ionosphere  parameters  (Ll  cycles) 


The  observables  are  replaced  by  its  residuals  as  shown  below. 
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The  geometry,  or  range  model,  and  the  ionosphere  model  are  replaced  by  its  variations  about 


the  a  priori  value. 
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If  an  a  priori  ionospheric  contribution  is  also  removed  from  the  prefit  phase  residuals,  then  the 
parameter  k  gets  replaced  by 


fc  =  Md,0n|0 


(4.25) 


The  next  step  is  to  further  expand  x2  and  to  form  the  normal  equations  for  least  squares 
analysis.  Reinserting  the  biases,  matrix  (  V  J  becomes 
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The  bias  parameter,  like  the  other  geometric  parameters,  can  be  rewritten  as 
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Now  [  V  ]  can  be  rewritten  with  explicit  paruals. 
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It  should  be  noted  that  the  parameters  (x,  y,  nu,  iil:)  above  represent  a  vector,  therefore,  the 
partial  derivatives  also  represent  vectors. 
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In  practice,  the  L 2  phase  ambiguity  parameters  are  not  directly  solved  for.  The  phase 
ambiguity  of  the  beat  frequency  of  LI  and  L2  signals  is  calculated.  The  beat  frequency  of 
approximately  350  MHz  is  obtained  by  differencing  the  two  carrier  frequencies.  The 
wavelength  of  this  beat  frequency  is  86  cm,  and  for  codeless  receivers  the  wavelength  is 
reduced  to  1  j  cm.  This  bias  known  as  the  wide  lane  bias.-*  It  is  determined  by  differencing  the 
LI  and  L2  bias  parameters.  In  liquation  4.29  the  n^  bias  is  replaced  by  the  wide  lane  bias 
(nL2  -  nti )•  This  requires  that  the  elements  of  the  transformation  matrix  in  Equation  4.29,  and 
only  this  equation,  be  replaced  by  the  following  expressions: 


1 1 1  —  1 1 1  +  1 1 2  (4.32) 

l:i  =  l2i  +  l22 .  (4.33) 


Equation  4.29  can  be  further  simplified  to: 
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or  in  a  more  compact  notation: 
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i  P  ]  -  partial  derivative  matnx 

z  -  parameter  vector 

Ot  -  transformed  observation  vector 


(4.35) 


.  (4.34) 


Reinserting  [  V  ]  back  into  Equation  4.1 1,  and  expanding  we  get: 
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Now,  x 2  can  be  easily  minimized  with  respect  to  the  adjustable  parameters  (  z  ).  This  is 
accomplished  by  taking  the  partial  derivative  of  Equation  4.36  with  respect  to  z,  and  setting  the 
equation  to  zero.  The  resulting  set  of  linear  equations  are  known  as  the  normal  equations,  and 
are  solved  for  z. 

[Ajz=[B]  (4.37) 
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Success  of  this  approach  depends  on  the  accuracy  of  the  constraint  placed  on  the  ionosphere. 
The  effect  of  the  ionospheric  contribution  are  inseparable  from  the  integer  ambiguity.  To  get 
accurate  integer  bias  values,  the  ionosphere  must  be  accurately  constrained.  If  the  ionospheric 
constraint  (  Oj0n  )  is  small,  the  ionospheric  parameters  cannot  be  sufficiently  adjusted  to  agree 
with  the  observations.  If  the  ionospheric  constraint  is  too  large,  the  uncertainty  is  passed  onto 
the  estimates  of  the  bias  parameters  and  it  becomes  difficult  to  fix  biases. 

4.2  Addition  of  the  C/A  code 

This  concept  of  estimating  ionospheric  parameters  can  be  extended  to  include 
pseudorange  observables,  in  particular  the  C/A  code.  Although  the  pseudorange  observable  is 
noisy  and  imprecise,  it  contains  unambiguous  range  information,  which  is  not  available  with 
phase  observable.  A  simple  example  can  show  the  advantages  of  adding  pseudorange 
information.  This  example  looks  at  the  effect  of  combining  pseudorr>m»e  and  phase 
observables,  with  the  ionosphere  effects  neglected.  Chi-squared  is  written  as: 
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°C/A  °L1  .  (4.42) 

Oc/a  -  C/A  code  observable  (LI  cycles) 

Gc/a  *  assumed  u  priori  standard  deviation  for 
the  C/A  code  observable  (LI  cycles) 


For  this  study  C/A  code  observations  have  been  convened  from  units  of  delay  to  phase. 
Minimizing  x2  with  respect  to  geometry  (dRCom)  and  bias  (nu)  parameters,  and  making  use  of 
Equations  4.23  and  4.27,  the  normal  equation  can  be  written  as: 


(4.43) 
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By  taking  the  inverse  of  (A) ,  we  get  a  covariance  matnx  of  errors  for  the  parameters. 

In  this  example  the  covariance  matrix  contains  elements  which  are  the  inverse  of  elements  in 
matrix  [A].  The  magnitude  of  elements  in  the  covariance  matrix  can  be  reduced  by  decreasing 
oLi  or  by  introducing  a  nonzero  Cc/A-  The  value  of  covariance  matrix  elements  can  also  be 
reduced  by  decreasing  Oc/a-  The  analysis  shows  that  for  a  given  geometry  (i.e.  the  pamals  are 
the  same),  the  C/A  code  or  other  pseudorange  information  will  help  reduce  the  uncertainty  in 
the  estimate  of  the  bias.  If  the  pseudorange  information  is  very  noisy  ( large  Oc/\  ),  the 
improvement  is  negligible. 

The  approach  used  in  Section  4. 1  can  expanded  to  include  an  additional  type  of 
observable.  In  this  study,  the  aim  is  to  add  the  ability  to  process  C/A  code  observables  along 
with  the  two  phase  observables.  First,  x~  must  be  reformulated  to  include  the  residual  of  the 
C/A  observable: 
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The  ionosphere  contribution  term  (k)  is  eliminated  as  previously  described.  The  difference 
now  is  that  there  are  three  observables,  and  we  need  to  transform  them  into  two  observables. 
The  transformation  becomes  a  fit  of  three  observable  to  the  desired  two  observables,  instead  of 
a  simple  rotation,  or  linear  combination,  of  just  two  observables.  To  account  for  any 
differences  in  X~  between  the  form  in  Equation  4.46  and  the  quadranc  form,  an  addinonal  term 

(c)  is  introduced  as  shown  below: 


Xj  =  [  V  ]  T[  W  ]  [  V  ]  +c 


(4,47) 


c  -  term  to  account  for  difference  in  x2  between 
Equations  4.46  and  4.47 

Now  the  transformed  observables  are  determined  as  follows: 
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Since  the  quantities  being  fitted  are  observables,  the  resulting  discrepancy  (c)  between  the  two 
forms  of  x2  will  again  only  be  a  function  of  the  observables.  None  of  the  parameters  will  be 
involved,  implying  that  when  y}  is  minimized  with  respect  to  the  parameters,  the  additional 
term  (c)  will  not  affect  the  normal  equations.  The  weight  and  transformation  matrices  are 
calculated  as  described  in  Section  4.1.  All  the  terms  are  reinserted  into  Equation  4.47,  which  is 
then  subtracted  from  Equation  4.46,  and  extra  term  (c)  can  now  be  determined.  As  before,  if 
the  wide  lane  biases  (n^  -  nu)  are  desired.  Equations  4.32  and  4.33  must  be  substituted  into 
Equation  4.34.  These  calculations  were  also  carried  out  by  Dr.  Gourevitch  using  Macsyma. 
Complete  expressions  for  the  weight  and  transformation  matrices,  and  the  term  (c)  is  given  in 
Appendix  B.  By  setting  the  weight  (  1/Cc/a2  )  for  the  C/A  code  to  zero,  it  is  apparent  that  the 
LI  and  L2  only  case  is  a  subset  of  this  formulation.  The  term  (c)  is  zero  if  the  weight  for  any 
observable  is  zero. 

The  normal  equations  are  formed  as  describe  in  Section  4  1  and  used  to  solv  e  for  the 
bias  parameters.  The  C/A  code  not  only  provides  range  information,  but  also  information  on 
the  ionospheric  contribunon.  The  sign  for  ionospheric  contribution  is  different  for  the  phase 
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and  pseudorange  observables  since  the  ionosphere  advances  the  phase  of  the  carrier  wave,  but 
retards  the  code  modulation. 

4.3  Simulation  of  C/A  code  observable 

To  study  the  effect  of  the  processing  C/A  code  and  phase  observations,  C/A  code 
observations  were  simulated  and  processed  with  real  phase  observation.  The  details  of 
forming  C/A  observations  are  discussed  in  this  section. 

Several  error  sources  corrupt  the  measurements  of  pseudorange  observations.  The  two 
major  sources  of  error  are  the  ionosphere  and  multipath.  The  C/A  code  observables  were 
created  by  modifying  real  LI  phase  observable  in  two  steps.  First,  the  position  of  one  station 
on  the  baseline  was  adjusted  using  the  LC  observable.  Since  this  observable  was  free  of 
ionospheric  contributions,  the  resulting  position  adjustment  was  also  free  of  ionospheric 
effects.  The  adjusted  position  was  then  used  to  calculate  the  theoretical  range  to  the  satellites. 
This  updated  theoretical  range  was  used  to  form  the  prefit  residuals,  which  were  used  to 
simulate  the  C/A  code  observations. 

The  ionosphere-free  position  adjustment  was  used  to  form  the  LI  and  L2  phase 
residuals.  These  residual  were  then  used  to  form  a  LC  observable,  as  shown  in  Equation  4.9. 
The  C/A  code  was  formed  as  shown  below: 

0C/A  =  2$^  -  0L1  +  (ionospheric  contribution  bias »  +  ( noise)  (4  49) 

=  ^Ll  “2  — ^  (  $L2~  8  ^L1  ) 

1  -  g 

+  (ionospheric  contribution  bias)  +  (noise)  ^  ^ 

The  ionospheric  contribution  \>.us  simulated  by  the  measuring  the  change  in  the  ionospheric 
contribution  relative  to  the  first  ep<vh,  using  LI  and  L2  phase  measurements  The  resulting 


30 


ionosphere  in  the  C/A  code  observable  has  the  opposite  sign  than  the  ionospheric  contribution 
in  the  phase  observables.  The  International  Reference  Ionosphere  (  IRI)  Model  was  used  to 
determine  the  mean  value  of  the  ionospheric  contribution.  This  model  is  capable  of  predicting 
spatial  and  temporal  variations  within  the  ionosphere.1  There  are  two  ways  in  which  the  model 
can  be  utilized  to  calculate  the  constant  term.  The  first  method  is  to  compute,  from  the  LRI 
model,  the  ionospheric  delay  at  the  first  epoch,  in  each  series  of  observations.  The  second 
method  involves  determining  the  constant  term  at  every  epoch,  by  subtracting  the  (L2  -  g  LI) 
term  from  the  value  determined  from  the  model,  and  than  calculating  the  average  value  of  this 
constant  over  the  entire  series.  The  first  method  was  used  in  this  study  to  determine  the 
constant  value. 

The  second  major  source  of  error  contaminating  the  C/A  code  observable  is  multipath. 
As  stated  in  Section  2.4.4,  multipath  propagation  errors  exhibit  a  period  of  approximately  15 
minutes.  This  period  is  slightly  more  than  twice  the  observation  interval  of  7  minutes.  This 
fact  allows  the  correlation  between  successive  observations  to  neglected,  and  permits  each 
observation  to  have  independent  statistics.  For  this  study  the  double  difference  C/A  code 
observations  were  assumed  to  have  uncorrelated  gaussian  noise  of  2  m  rms.  The  one-way  C/A 
code  observations  were  modeled  with  1  m  of  random  gaussian  noise.  In  forming  differences, 
the  variances  of  one-way  noise  terms  add,  increasing  the  noise  present  on  double  difference 
observables. 

The  receiver  system  noise  (kT)  was  considered  negligible.  Tropospheric  contribution 
was  also  considered  negligible  for  this  study. 

The  value  of  chi-square  per  degree  of  freedom,  or  reduced  chi-square  (  Xv:h  was  used 
to  adjust  the  assumed  a  priori  measurement  errors  for  the  phase  observables  and  the 


1  Rjwer,  K  ,  D  Bilitza,  Study  of  Ionospheric  mil  Tropospheric  Models",  S'ASA  Report  V.Sft  1V85. 
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ionosphere,  to  assure  that  they  are  properly  weighted  relaove  to  the  C/A  code.  Chi-square  can 
also  be  written  in  terms  of  estimated  and  assumed  a  priori  measurement  errors. 1 
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s  -  estimated  measurement  error 

__  *> 

a  "  -  average  assumed  a  pnori  measurement  error 

v  -  number  of  degrees  of  freedom 


The  reduced  chi-square  (  Xv2 )  can  be  defined  as: 
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(4.54) 


Xv2  ■  chi-square  per  degree  of  freedom 


If  the  estimated  measurement  errors  accurately  represent  the  a  prion  measurement  errors,  the 
value  Xv2  is  approximately  1 .  If  the  a  prion  errors  resulted  in  Xv2  not  equal  to  one,  these  a 
priori  errors  (a, )  were  scaled. 


^  Bcvington.  P  R  ,  Dmq  RedtM  lion  and  f  rrt,r  A  ruivtif  r  O u*  Rh  » ' .1  dl  Si  wn*  ei .  SU1  »r  H.  II  Book  Co  JSc  * 
York.  NY.  pp.  187-189.  1%9 


(4.55) 

-  initial  assumed  a  priori  measurement  error 

-  scaled  assumed  a  priori  measurement  error 

forcing  Xv2  to  equal  one.  The  a  priori  error  for  the  ionosphere  was  also  scaled  to  preserve  the 
relative  weighting  of  the  phase  observables  and  the  ionosphere.  Since  the  error  in  the  C/A  code 
was  simulated,  its  a  prion  measurement  error  was  accurately  known.  Phase  and  ionospheric 
contribution  a  priori  errors  were  scaled  to  assure  that  the  relative  weighting  of  the  phase 
observables,  C/A  code  observables,  and  the  ionospheric  delay  were  consistent. 

LI  and  L2  phase  observations,  recorded  November  8,  1988  in  east  Texas,  were  used  to 
form  C/A  code  observations.  Observations  for  baseline  lengths  ranging  from  10  km  to  330  km 
were  recorded.  For  this  study,  observations  from  10  km,  100  km,  and  330  km  baselines  were 
processed.  The  initial  a  prion  measurement  and  ionosphere  uncertainties  were  based  on 
previous  analysis  of  this  data.1  The  observations  were  processed  for  single  baselines.  The 
simulated  C/A  code  observations  were  processed  with  the  phase  observables,  having  scaled  a 
prion  errors.  To  keep  the  analysis  simple,  only  one  station  position  and  the  LI  and  L2-L1  bias 
parameters  were  solved  for.  These  results  were  compared  to  solutions  performed  without  the 
C/A  code. 
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Carrier  Phase  Ambiguity  ”.  Presented  at  the  f-ifth  Injfncdioncil  Geodetic  Symposium  on  Satellite  Sc* 

Mexico  State  University,  Las  Cruses.  March  1  '  P.  PKd 


33 


5.0  Results 


This  section  presents  test  results  which  establish  the  validity  of  the  compute  •  programs, 
and  show  the  effect  of  combining  C/A  code  observations  with  LI  and  L2  phase  observations. 
First  the  validity  of  the  program  was  established,  using  simulated  LI  phase,  L2  phase,  and  C/A 
code  observations.  After  verification  of  the  programs,  simulated  C/A  code  observations  were 
processed  with  real  phase  observations. 

5.1  Verification  of  programs 

Upon  completion  of  modifications  to  the  software,  the  computer  programs  were  tested 
using  two  separate  procedures.  No  ionospheric  effect  was  modeled  for  the  simulated 
observations.  The  a  priori  ionospheric  parameter  (k)  assumed  to  have  a  zero  mean.  It  was  not 
fixed  and  allowed  to  vary  wtith  with  a  standard  deviation  of  0.001  ppm.  A  set  of  simulated 
observations  were  created,  and  processed  through  the  least  squares  algorithm.  In  the  simulated 
observations,  noise  was  added  by  using  gaussian  random  number  generator,  fhe  nnse  for  the 
double  difference  phase  observations  had  a  zero  mean,  with  a  standard  deviation  of 
approximately  0.2  LI  cycle  of  phase,  or  an  equivalent  path  length  of  2  cm.  The  noise  for  the 
double  difference  C/A  code  observations  also  had  a  zero  mean,  with  a  standard  deviation  of 
approximately  20  LI  cycles  of  phase,  or  an  equivalent  path  length  of  2  m.  The  prefit  residual 
for  each  observable  was  simulated.  The  prefit  residuals  were  simulated  for  two  baselines:  a  10 
km  baseline  and  a  100  km  baseline.  The  simulations  were  bared  on  a  actual  observation 
schedule.  In  processing  the  simulated  observations,  13  parameters  were  adjusted:  3  station 
coordinates,  5  LI  bias  parameters,  and  5  L2-L1  (wide  lane)  bias  parameters. 

The  equations  were  initially  checked  for  consistency.  The  assumed  a  priori  standard 
deviation  of  the  observations  must  be  con  istent  with  the  standard  deviation  of  the  simulated 
random  noise.  This  check  was  accomplished  by  looking  at  value  of  Chi-square  per  degree  of 
freedom  (  %v2  )•  If  the  equations  are  correct,  and  consistent  weights  are  used,  ihe  value  for 
Chi-square  per  degree  of  freedom  should  he  approximately  one. 


Since  the  simulated  errors  of  observations  were  noise  with  zero  mean,  the  adjustments 
to  the  parameters  should  also  be  random  variables  with  zero  mean.  The  uncertainty  of  the 
adjusted  parameters  will  also  be  random  variables.  Several  factors  influence  the  parameter 
adjustments.  Averaging  of  many  observations  reduces  the  size  of  adjustments.  The  geometry 
of  the  satellite  spacing  in  the  sky,  tends  to  increase  the  parameter  adjustments,  since  in  general 
the  satellites  are  not  evenly  distributed  across  the  sky. 

For  each  baseline  the  observations  were  processed  in  two  groups.  First  the  LI  and  L2 
phase  residuals  were  processed,  and  then  C/A  code  residuals  were  added  and  processed.  The 
results  for  the  10  km  baseline  are  given  in  Table  5.1 .  The  uncertainties  given  below  are 
defined  as  three  times  the  formal  standard  deviation  of  error  of  the  related  parameters  and  is 
multiplied  by  yv2. 

The  adjustment  of  the  position  coordinates  and  the  LI  bias  parameters  are  zero  to 
w  ithin  a  few  millimeters.  The  uncertainty  for  the  station  coordinates  and  the  LI  bias 
parameters  are  on  the  order  of  1  cm,  which  was  the  assumed  a  priori  standard  deviation  of  the 
observations.  The  uncertainty  in  the  L2-LI  bias  is  on  the  order  of  5  cm.  The  L2-L1  bias 
uncertainty  is  greater  than  that  for  the  LI  bias,  since  it  is  a  linear  combination  of  the 
uncertainties  in  the  LI  and  L2  bias.  It  was  assumed  that  the  observations  were  from  a  codeless 
receiver.  The  wavelength  of  the  simulated  LI  cycle  phase  corresponds  to  1/2  LI  phase,  and 
the  simulated  L2-L1  cycle  phase  corresponds  to  1/2  (L2-L!)  phase.  A  phase  change  of  1  cycle 
at  LI  frequency  and  L2-L1  beat  frequency  is  equivalent  to  a  change  in  path  length  of  9.5  cm 
and  43  cm  respectively.  Similar  results  were  obtained  for  the  100  km  baseline,  and  are  not 
presented  here. 


Table  5.1  10  km  baseline  parameter  adjustments  and  uncertainties 


v  2 

Parameters 

LI  +  L2 

LI  + 

L2  +  C/A 

1 

Adjustment 

0.926 

Uncertainty 
(3  *CTUnceruimy*Xv't) 

Adjustment 

0.903 

Uncertainty 

(  1  *  y  * 

■ J  ^  uncertainty  Av 

Latitude  (rad) 

0.143xl0'7 

0.49x10" 7  (0.005m) 

0.143xl0"7 

0.53x10" 7  (0.005m) 

Longitude  (rad) 

0.186xl0"7 

O.HxlO"6  (0.013m) 

0.189xl0"7 

0.53X10"6  (0.013m) 

Radius  (km) 

-0.1 16xl0'5 

0.12x10"*  (0.012m) 

-0  113xlO"5 

0. 53x10"*  (0.012m) 

LI  Bias  (LI  cycles  ): 

1 

-0.0439 

0.101 

-0.0443 

0.0995 

i 

-0.0322 

0  115 

-0.0327 

0.114 

3 

-0.0587 

0.113 

-0.0588 

0.112 

4 

-0.0032 

0.120 

-0.0041 

0.118 

5 

-0.0611 

0.169 

-0.0594 

0.167 

L2-L1  Bias(L2-Ll  cycles): 

1 

0.0531 

0.120 

0  n.538 

0.119 

2 

0.0240 

0.126 

0(i248 

0.124 

3 

0.0802 

0. 1 30 

0  UX07 

0.128 

4 

0.0587 

0  1 19 

0  0603 

0.117 

5 

0.0849 

0.138 

0.0824 

0.137 

5.2  Integer  ambiguity  determination 

Next,  the  effect  of  including  C/A  code  observables  on  the  uncertainties  in  the  estimation 
of  station  position  coordinates  and  the  bias  parameters  was  studied.  Observations  for  a  given 
baseline  were  processed  with  and  without  the  simulated  C/A  code  observations.  The  a  prion 
standard  deviation  of  double  difference  phase  observauons  and  the  ionosphere  are  listed  in 
Table  5.2.  The  standard  deviation  was  calculated  by  taking  the  constant  term  and  adding  a 
second  term,  that  is  proportional  to  baseline  length. 


Tabic  5.2  Assumed  a  priori  standard  deviation  of  double  difference  observations 

Constant  term  Distance  proportional  term 

Phase  observations  3.5  mm  1.0  ppm 

Ionospheric  contribution  0.0  mm  30.0  ppm 

Observations  for  10  km,  100  km,  and  330  km  baseiines  were  processed.  On  the  10  km 
baseline,  the  addition  of  C/A  code  observations  resulted  in  a  slight  increase  in  the  uncertainty  of 
the  parameter  estimate,  as  shown  in  Table  5.3.  In  both  cases  the  uncertainty  of  the  bias 
parameters  were  small  enough  to  enable  these  parameters  to  be  fixed  to  an  integer  value,  using 
the  method  described  in  Footnote  1  given  below. 

One  possible  explanation  for  results  shown  in  Table  5.3  is  the  erroneous  assumption  of 
zero  a  priori  ionospheric  contribution  on  both  the  phase  and  pseudorange  observables.  This 
assumption  was  implicitly  made  by  setting  dion  to  zero  (see  Section  4.1).  The  ionospheric 
effect  (k)  on  LI  double  difference  observable  was  also  assumed  to  be  an  independent  random 
variable  with  zero  mean  and  a  standard  deviation  (chon)  which  was  a  function  of  baseline 
length.  In  reality  the  observations  contain  a  non-zero  a  prion  ionospheric  contribution,  which 
is  correlated  in  time.  This  systematic  error  was  also  present  in  the  C/A  code;  therefore,  by 
adding  pseudorange  observations  to  phase  observations  the  parameter  estimates  are  further 
corrupted. 


Table  5.3  10  km  baseline  uncertainties  in  parameter  estimates 


Parameters 

LI  +L2 

LI  +  L2  +  C/A 

Uncertainty 

0*Ouncerumty*Xv7) 

Uncertainty 

(  3  *  O  uncertain  ty  *  Xv  7  ) 

Latitude  (rad) 

0.31x10- 7  (0.003m) 

0.32xI0'7  (0.004m) 

Longitude  (rad) 

0.86x10' 7  (0.008m) 

0.88xl0'7  (0.009m) 

Radius  (km) 

0.78xl0‘5  (0.008m) 

0.80xl0‘5  (0.008m) 

LI  Bias  (LI  cycles  ): 

1 

0.407 

0.414 

2 

0.425 

0.433 

3 

0.439 

0.447 

4 

0.401 

0.408 

5 

0.465 

0.474 

L2-L1  Bias  (L2 

■LI  cycles): 

1 

0.117 

0  119 

2 

0.122 

0.124 

3 

0.126 

0.128 

4 

0.115 

0.117 

5 

0.132 

0.135 

Results  from  the  100  km  baseLine  are  shown  in  Table  5.4.  The  addiuon  of  C/A  code 
observations  on  the  100  km  baseline  caused  a  slight  decrease  in  the  uncertainty  of  the  parameter 
estimates,  compared  with  the  solution  using  only  phase  observations.  For  the  10  km  baseline, 
the  parameter  estimate  uncertainties  did  not  decrease.  The  large  uncertainties  of  the  bias 
parameters  prevented  any  biases  from  being  fixed  to  an  integer  value. 


Table  5.4  100  km  baseline  uncertainties  in  parameter  estimates 


Parameters 

LI  +L2 

LI  +L2  +  C/A 

Uncertainty 

(  3  *  ^uncertainty  *  Xv^ ) 

Uncertainty 

( 3  *  CuncctLauity  *  Xv^  ) 

Latitude  (rad) 

O.BxlO-6  (0.015m) 

0.13x1 0"^  (0.015m) 

Longitude  (rad) 

O  ^rlO-6  (0.035m) 

0.35X10"6  (0.034m) 

Radius  (km) 

0.32x10^  (0.032m) 

0.31x10^  (0.031m) 

LI  Bias  (LI  cycles ): 

1 

1.75 

1.72 

2 

1.81 

1.78 

3 

1.85 

1.82 

4 

1.75 

1.72 

5 

1.94 

1.90 

L2-L1  Bias(L2-Ll  cycles): 

1 

0.502 

0.492 

2 

0.519 

0.508 

3 

0.530 

0.520 

4 

0.501 

0.491 

5 

0.552 

0.541 

Results  from  the  330  km  baseline  were  qualitatively  similar  to  the  10  km  baseline 
results,  but  the  magnitude  of  the  uncertainties  themselves  were  much  larger.  The  results  are 
not  presented  here.  The  coordinate  uncertainties  were  several  centimeters,  and  the  bias 
parameter  uncertainties  were  approximately  4  cycles  for  the  L2-L1  biases  and  around  14  cycles 
for  the  LI  biases.  No  biases  were  able  to  be  fixed  for  this  baseline. 

To  account  for  the  lack  of  improvement  from  the  addition  of  C/A  code  observables, 
further  tests  were  performed  to  examine  the  accuracy  of  the  C/A  code  simulation.  The  standard 
deviation  of  random  noise  for  the  simulated  double  difference  C/A  code  observ  ations  was 
reduced  from  2  m  to  20  cm,  and  the  observations  were  processed  for  the  10  km  baseline. 
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The  resulting  solution  showed  an  increase  in  the  uncertainties,  instead  of  a  decrease. 

The  value  of  xv2  increased  from  approximately  I ,  with  2  m  rms  noise,  to  3.27.  This  increase 
in  Xv2  indicates  that  the  assumed  20  cm  assumed  a  priori  standard  deviation  for  the  CJA  code 
was  less  than  the  actual  standard  deviation  of  the  simulated  observations.  By  reducing  the 
assumed  a  priori  standard  deviation  on  the  C/A  code,  the  weight  of  these  observations  were 
greater,  allowing  the  unmodeled  sources  to  corrupt  the  baseline  measurement.  Using  C/A  code 
observations  with  2  m  rms  noise,  and  setting  the  assumed  a  priori  standard  deviation 
accordingly,  significantly  reduces  the  weight  of  the  pseudorange  observations  relative  to  the 
phase  observations.  The  lower  weight  reduces  the  effect  of  the  unmodeled  errors  in  the  C/A 
code  data  on  the  baseline  measurement. 

The  role  of  the  ionospheric  contribution  was  further  studied,  to  account  for  results 
presented  above.  Although  the  effect  of  the  ionosphere  on  pseudorange  and  phase  observables 
is  equal  and  opposite,  the  ability  to  retrieve  information  on  the  ionosphere  from  each 
observable  differs.  Due  to  integer  ambiguity  present  in  the  double  difference  phase  observable, 
it  is  difficult  to  get  a  accurate  estimate  for  the  ionospheric  parameter,  without  fixing  the  biases. 
Pseudorange  observables  do  not  contain  an  ambiguity  parameter,  making  it  possible  to  get  an 
accurate  estimate  of  ionospheric  parameter  (k).  The  C/A  code  observations  are  only  available  at 
a  single  frequency,  reducing  its  potential  rvr  estimating  the  ionospheric  parameter. 

Pseudorange  observations  at  both  the  LI  and  L2  frequencies  would  allow  a  more  accurate 
estimation  of  the  ionospheric  contribution. 

Further  test  cases  were  run  to  determine  the  viability  using  the  ionospheric  information 
on  the  C/A  code  to  constrain  the  estimation  of  the  ionosphere.  This  was  tested  by  loosening 
the  constraint  on  the  ionosphere  (i.e.  increasing  oIon  )  to  1000  ppm,  on  the  10  km  baseline. 
Phase  observations  were  processed  with  pseudorange  observations,  having  2  m  and  20  cm  of 
rms  random  noise.  The  two  test  cases  are  listed  in  Table  5.5. 
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Table  5.5  Assumed  a  priori  standard  deviation  for  two  test  cases 


Case  1  Case  2 


Phase  observations  2.9  mm 

C/A  code  observations  2000  mm 

Ionospheric  contribution  1000  ppm 


2.9  mm 
200  mm 
1000  ppm 


Using  C/A  code  observations  with  2  m  rms  random  noise,  parameter  estimates  were  similar  to 
those  obtained  using  only  LI  and  L2  phase  observations  with  a  tighter  ionospheric  constraint 
(8  ppm).  This  can  be  attributed  to  the  small  weight  given  the  C/A  code  relative  to  phase 
observations.  The  value  of  Xv2  was  1.04.  The  uncertainty  of  the  integer  bias  parameters 
increased  by  a  factor  of  100,  reflecting  the  large  a  priori  ionospheric  contribution  uncertainty. 
The  large  a  priori  uncertainty  implies  that  the  LI  and  L2  phase  observables  are  being  combined 
into  an  LC  observable,  resulting  in  bias  parameters  which  are  highly  correlated.  For  the 
second  case,  C/A  code  observations  with  20  cm  of  measurement  error  were  processed  with 
phase  observations.  It  was  hoped  that  in  this  case,  by  increasing  the  weight  of  the 
pseudorange  observable,  the  ionospheric  error  estimate  would  improve.  Instead,  results 
similar  to  previous  case  were  obtained,  except  that  the  uncertainties  in  the  position  coordinates 
were  larger  by  50  %  and  Xv2  increased  to  3.25.  These  results  indicate  that  the  C/A  code 
observable  does  not  improve  estimation  of  the  ionospheric  error. 
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6.0  Conclusions 

This  study  indicates  that  addition  of  the  C/A  code  does  not  improve  the  estimation  of 
baseline  parameters.  This  statement  must  be  taken  cautiously,  since  there  are  several  sources 
of  error  in  this  study.  Two  primary  sources  of  weakness  are  the  lack  of  an  a  priori  model  for 
the  ionosphere,  and  the  simulation  of  C/A  code  observations.  The  lack  of  an  a  priori , 
parameterized,  ionosphere  model  prevented  the  ionospheric  contribution  from  being  more 
thoroughly  removed  from  the  observations.  The  use  of  an  imperfect  ionosphere  model  for  the 
mean  value  of  the  ionosphere,  introduced  errors  in  the  simulated  C/A  code,  which  were 
inconsistent  with  the  phase  observations.  It  has  proved  to  be  difficult  to  make  the  simulated 
C/A  code  observations  as  realistic  and  consistent  with  the  phase  ob^rvations. 

It  appears  that  pseudorange  information  is-redundant  for  short  baselines,  since  the 
integer  ambiguity  parameters  can  determined  for  such  baselines.  The  pseudorange  information 
is  more  valuable  on  the  longer  baselines,  where  it  is  difficult  to  fix  the  integer  ambiguity 
parameters. 
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Appendix  A 


below. 


The  elements  of  the  weight  matrix,  for  LI  and  L2  phase  only  algorithm,  are  given 


[  W  ]  = 


wn  w  j  2 
W  2  1  w22 


W|  = 


1 

w3  =  — 


(A.l) 


(A. 2) 


(A. 3) 


(A. 4) 


(A. 5) 


4  2 

w2(w1  +  w4)g  +W;(w4-2  w2)g  +  w,w2 

wn  _  - 

(Wl  +  w4)g  +w2  (A6) 

W4  (W|  +  w:  )  g“ 

W|2  = - ; - 

fw,  +  w4)g‘  + w2  (A  .7) 

w12  =  w  2 1  (A. 8) 
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w22  = 


W4  (w  1  g2  +  W2) 
(W,  +  w4)g2+  w2 


(A. 9) 


The  elements  of  the  transformation  matrix,  for  LI  and  L2  phase  only  algorithm,  are  given 
below. 


[T]  = 


'll  l12 

'2 1  l22 


=  - 2 


1-g 


2  =  '  2  2 


1  -  g 


t21  =  - 


1-g 


'22  = 


_  g 


-  g 


(A.  10) 


(A.ll) 


(A. 12) 


(A.  13) 


(A.  14) 
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Appendix  B 

The  elements  of  the  weight  matrix,  for  LI  phase,  L2  phase,  and  C/A  code  algorithm, 
are  given  below. 


[  w  ]  = 


w,, 

w  1  2 

w  2 1 

w  22 

(B.l) 


The  the  weight  terms  (wIt  w2,  W3,  and  w4)  are  defined  in  Appendix  A,  Equations  A. 2  through 
A. 5. 


w 


4  -,2 

W2(w,+W3+W4)  g  +  (w,+  w3  )(2w2+w4  )+4w,w3  g+w2(w!+w3) 

n  -  ; - J - 


(W,  +  W3+W4  )g  +  w 2 


(B.2) 


w12 


w4  (w  1  -f  w2  •  w3)  gz 

(Wi  +  W3  +  W4)g2  +  W2 


W12  =  w  2  j 


W4  (W|+  W3  )  g‘  +  W;  W4 
(wi+W3  +  w4)g2  +  w2 


(B.3) 


(B.4) 


(B.5) 


The  elements  of  the  transformation  matrix,  for  LI  and  L2  phase  only  algorithm,  are  given 
below. 


Ml  M  2  l32 
( 2  j  t  2  ^2 


(B.6) 


Mi  = 


W|llWj-Wi)g'  +  W|Wj 


4 

w2(wj  +  w3)g  +  2  Wi(wj  -  Wj  )  +  2  w,  w3 


+  w 2  ( w  j  +  w 3  )  (B  ?) 
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t12  - 


w2(w,  +  w3)g  +  w  2  ( w ,  -  w  3)  g 


W2(w,+W3)g  +  2 


w2(w3-w1)  +  2w1w 


g%w,(w1+w3) 


*13  = 


w3  (  2  w,  +  W2  )  g  +  w2  w3 


W2  (w,  +  w3 )  g  +2 


vv2(w3-w1)  +  2w;w 


g  +  W  2(w1  +  w- 


*21  “ 


wj  w2g4  +  2  (w3-  w2  )  g2 


w2(wj  +  W3  )  g  +2 


Wj(w3  -  W|  )  +  2  Wj  w3 


g  +  w2(w,  +  w3  ) 


t22  - 


w2(2w3-w,)g  +w2(w,  +  w3)g 


W2  (W|  +  w3 )  g  +2 


w2(w3- w,  )  +  2  w,  w- 


g  +  w2(w,  +  w3) 


*23  ~  ' 


w2  w3  g  +  (  2  w,  +  w2)  g‘ 


W2(W(  +  W3)g  +2  W2(w3-w1)  +  2wtw3jg‘  +  w2(w1  +  w3) 


The  extra  term  (c)  in  Equation  4.47  is: 


W|  w  2  w  3 

0 

Ll(  1  +  g2)-2d)L23  -  0r/A  (  1  -g*) 

L  . 

w2(w,  +  w3  )  g4  +  2 

w2(w3-W|)  +  2w|W3 

2 

g  +  w  ;  (  W  j  +  W  3  ) 

Note  that  all  the  elements  of  the  transformation  matrix  and  the  extra  term  (c)  have  the 
denominator.  The  extra  term  is  zero  if  any  of  the  weights  are  equal  to  zero. 


( B  S ) 


(B.9) 


(B.  10) 


(B.ll) 


fB. 1 2) 


( B .  1 3 ) 

same 
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